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Abstract 
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Wc consider the tricritical Ising model on a strip or cylinder under the integrable perturbation 
by the thermal 1^91.3 boundary field. This perturbation induces five distinct renormalization 

■ group (RG) flows between Cardy type boundary conditions labelled by the Kac labels (r, s). We 
^ I study these boundary RG flows in detail for all excitations. Exact Thermodynamic Bethe Ansatz 
Q . (TBA) equations are derived using the lattice approach by considering the continuum scaling 

I limit of the A4 lattice model with integrable boundary conditions. Fixing the bulk weights to 

. their critical values, the integrable boundary weights admit a thermodynamic boundary field 

I ^ which induces the flow and, in the continuum scaling limit, plays the role of the perturbing 

. boundary field <pi,3- The excitations are completely classified, in terms of string content, by 

' (m, n) systems and quantum numbers but the string content changes by either two or three 

^ . well-defined mechanisms along the flow. We identify these mechanisms and obtain the induced 

10 ' maps between the relevant finitized Virasoro characters. We also solve the TBA equations 

. numerically to determine the boundary flows for the leading excitations. 

00 ■ 

^ ■ 1 Introduction 

o : 

' Quantum Field Theories with a boundary have received a lot of attention recently due to their 

■ applications in Condensed Matter, Solid State Physics and String Theory (D-branes). A problem 
Qh! of great interest is the Renormalization Group (RG) flow between different boundary fixed points 

of a Conformal Field Theory (CFT) that remains conformal in the bulk. Many interesting results 
have been achieved, and flows have been studied for minimal models and for c = 1 CFT (see e.g. 
• 1— I ! ^ and references therein). Numerical scaling functions for the flow of states interpolating two 
^ I different boundary conditions can be systematically explored by use of the approximate Truncated 
'• Conformal Space Approach (TCSA) EiEllll. 

A beta function can be defined for the boundary deformations, much the same as for the bulk 
perturbations of conformal field theories [S]. The conformal boundary conditions play the role of 
ultraviolet (UV) and infrared (IR) points of the flow. One flows away from the UV fixed point 
by perturbing with a relevant boundary operator and flows into an IR fixed point attracted by 
irrelevant boundary operators. 

Among the possible boundary perturbations of a CFT there are some integrable perturbations 
that preserve an infinite number of conservation laws. In this case the flows are amenable to 
investigation by exact methods such as commuting transfer matrices and Bethe ansatz techniques. 
One of the most celebrated of these methods is the Thermodynamic Bethe Ansatz (TBA) [H] giving 
a set of non-linear coupled integral equations governing the scaling functions along the RG flow. 
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Boundary TBA equations were first obtained with the usual scattering approach [71 0] for the 
groundstate, a few excited states and the boundary entropy in the minimal models. In particular, 
the groundstate and boundary eutopics of the Tricritical Ising Model have been studied [S] within 
this approach. 

In this paper we use a lattice approach to obtain exact TBA equations for all the excitations 
of the integrable boundary flows of the Tricritical Ising Model (TIM) with central charge c = jq. 
This model has interesting applications in Solid State Physics and Statistical Physics. Its Kac table 
of conformal weights h^^s is 
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We have drawn only half of the table, by taking into account the well known Z2 Kac table symmetry 
(r, s) = (4 — r, 5 — s). The number of independent chiral primary fields is thus six. 

Cardy-type conformal boundary conditions for minimal models with diagonal modular invari- 
ant partition functions of type {Ap,Ap') are in one-to-one correspondence with the chiral primary 
fields in the Kac table. So the TIM (A3, ^4) admits six different conformal boundary conditions 
which we denote by -B(r.,s)) = 2, 3 and s = 1, 2. 

Let us consider the TIM defined on a strip of width L in the two-dimensional (x, t) plane so 
that < X < L and —oo<t< +00. Let us impose boundary conditions and B(^^j^gi-^ on the 

left {x = 0) and right (x = L) edge respectively. We denote this situation by B(^j.,s)\(r',s')- Iii the 
sequel we are interested in boundaries of type (r, s)|(l, 1) which we denote by = B(^j.,s)\(i,i) = 
^{r,i)\{i,s)- For these boundary conditions the partition function reduces to a single character |^ 

^%,,)(9) = %s)|(l,l)(9) = Xr,s{q) (1-1) 

where Xr,siQ) denotes the character of the irreducible Virasoro representation labelled by (r, s) at 
central charge c = jq. 

For each given boundary condition ;B(r,s) there is a set of boundary operators ipu^v that live on 
the edge. If we want to keep the boundary condition unchanged all along the edge, these operators 
must be restricted to the conformal families appearing in the OPE fusion of the Virasoro family (r, s) 
with itself: {u,v) G {r,s) x (r, s). They are distinguished in terms of their conformal dimensions 
hr^s as relevant (/ir,s < 1) and irrelevant {hr^s > !)• Of course only relevant perturbations break 
scale invariance at the boundary in such a way to get out of the fixed boundary point of a specific 
boundary condition and flow to another one. So if we want to consider possible relevant boundary 
perturbations of the TIM, i.e. QFTs described by the action 

/+00 
dtipuA^ = 0,t) (1.2) 
-00 

where S(^r,s) denotes the action of TIM with boundary condition ^(^.^s), we have to restrict to the 
possibilities as shown in Tabled 

Actually there are two physically different flows for each "pure" perturbation (i.e. containing 
only one operator ipi^2 or ipi^s), flowing to two possibly different IR destinies. This is achieved 
by taking different signs in the coupling constant in the case of 991^3 flows, and real or purely 
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Tabic 1. 



boundary condition 


boundary perturbations 


^(1,1) = '^(3,4) 




none 


^(2,1) = ■^(2,4) 





none 


^(3,1) = ■^(1,4) 


+ 


none 


^(1,2) = -6(3,3) 


-0 


fl,3 


^(2,2) = •6(2,3) 


d 


9^1,2, ^1,3 


^(3,2) = ■6(1,3) 


0+ 


•^1,3 



imaginary coupling constant in case of (fi^2 ones. The boundary condition ^(2,2) can be perturbed 
by any linear combination of the fields (^1,2 and V9i,3. The symbols +, — , 0+, —0, 0, d represent other 
ways lllj to denote the TIM conformal boundary conditions and we give them only to make 
contact with the existing literature. 

Along the physical flow from UV to IR, the boundary entropy associated to each ^(r,s) 

py/^ sin f sin f 
\5J ^sm f sm f 

decreases [S], so we only expect flows between boundary conditions where the initial boundary 
entropy is larger than the final boundary entropy. The possible conformal boundary conditions have 
been studied by Chim ^0] and the flows connecting them by Affleck The picture is summarized 
in Figure n Integrability can be investigated in a manner similar to the bulk perturbations and it 
turns out that the flows generated by pure ipi^s and ipi^2 perturbations are integrable. In contrast, 
the flow starting at -6(2.2) a perturbation which is a linear combination of (pi^2 and 931,3 is 
strongly suspected to be non-integrable. Notice the Z2 symmetry of Figure which is related to 
the supersymmetry of the TIM. Investigation of the supersymmetric aspects of the TIM boundary 
flows is, however, beyond the scope of the present paper. 

In this paper we derive exact TBA equations for all excitations of the integrable boundary flows 
of the TIM from the functional relations for the transfer matrix of the lattice RSOS model |12j 
which is a member of the Al Andrews-Baxter- Forrester models JHl- Using general techniques jl41 
I15j of the lattice approach for turning functional equations into non-linear integral equations, we 
derive the TBA equations for all excitations and solve them numerically to determine the flows 
for the leading excitations. Our preliminary results for the flow xi,2 ^ X2,i were announced in 
|16j . Although we only consider the TIM here, our approach is quite general and should apply, for 
example, to the integrable boundary flows of all minimal models. A similar lattice approach has 
been successfully applied ^Zj to the massive and massless bulk thermal KG flows of the TIM. We 
stress that we are not addressing here the interesting question of predicting the allowed pattern of 
flows between fixed points. Rather, given a particular flow, the thrust of the lattice approach is to 
calculate the flow of all excitations and to predict, excitation by excitation, where these levels flow. 

The layout of this paper is as follows. In Section 2, we discuss the TIM as the continuum 
scaling limit of the lattice model. We describe the integrable boundary conditions ^Sj that 
give rise, in the continuum scaling limit, to the conformal boundary conditions labelled by the Kac 
labels {r,s), the associated double row transfer matrices JH] and their functional equations. We 
summarize the results of O'Brien, Pearce and Warnaar [20] which apply at the conformal fixed 
point of the TIM before the boundary perturbation ^ = Im(^iatt) is turned on. In particular, we 
review the classification of the allowed patterns of zeros of the transfer matrix eigenvalues in terms 
of (m, n) systems and quantum numbers and summarize the TBA equations and their solutions for 
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Fig. 1. The boundary flows between TIM conformal boundary conditions. Pure (pi,^ and (pi^2 flows are 
integrable. More information can be found in The boundary entropy associated with the +&— boundary 
condition is Uncar in the Cardy type boundary entropies: (7(^fo_) — (7(i,i) +5(3,1). 



the finite-size energy spectra in terms of finitized characters. In Section 3, we consider the three 
boundary flows with constant r, namely, xi,2 > Xi,!; X3,2 ^ X3,i aiid X2,2 > X2,i- We describe 
two mechanisms A, B under which the patterns of zeros change during these flows and in each case 
we show that this is consistent with an explicit map between finitized characters. In each case, 
we also solve the functional equations to obtain the boundary TBA equations describing the flow. 
In Section 4, we repeat these considerations for the two boundary flows with variable r, namely, 
Xi,3 ^ X2,i and xi,2 ^ X2,i- In these cases we find three mechanisms A, B, C. In Section 5, we 
discuss the numerical solution of the boundary TBA equations for all five thermal boundary flows. 
Finally, we finish with some conclusions in Section 6. 

2 The tricritical Ising model 

The conformal unitary minimal model with central charge c = 7/10 describes the critical behaviour 
of the tricritical Ising model. A lattice realization of this universality class is given by the critical 
^4 ABF model on the critical line separating regimes III and IV. Since we are interested in working 
in the presence of a boundary we define, following [20], a double row transfer matrix with (ri,ai) 
and (r2,a2) boundary conditions by the following diagrammatic representation 



ri ri o'2 ""jv-i r2 r2 




ri ri CTi (T2 o-jv-i r2 r2 
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The bulk weights are fixed to their critical values 
W 











c 






it 






u 




b 


a 





sin(A — ti) sinu SaSc 

sm A sm A V Jb'Jd 



(2.2) 



where the heights a,b,c,d G {1,2,3,4} on any edge must respect the A4 adjacency rule. The 
physical range of the spectral parameter is < u < A where A = ^ is the crossing parameter. The 
crossing factors are Sa = sin a A/ sin A. 

The most general integrable boundaries for critical bulk weights are given by (3.28) in jl8j . 
They are labelled by a pair (r, a) and depend upon boundary spins b, c, d, the spectral parameter 
u and a boundary thermodynamic field ^ (interaction) 



d 



u, C 



(2.3) 



We only need to consider diagonal boundary weights in the sense that b = d. Explicitly, the relevant 
boundary weights for our purposes are 

c±l 



C± 1 
C± 1 



n, C 




c±l 



S(rTc+a)/2 S{c±aTr)/2 sin(C + u) Sm{Xr + £, - u) 

+5(,.±c+a)/2 5'(c=Fa±r)/2 sin(^ - u) sm{Xr + ^ + u) 



(2.4) 



sinA Sr (5c5c±i)^/^ cosh2Im(e) 

where r = 1, . . . , 3 and a = s = 1, . . . , 4. Here a and 6 = d = c it 1 must be adjacent at fusion level 
r, that is, a — b = ±1 for r = l, o + 6 = 4, 6 for r = 2 and a + b = 5 for r = 3. The normalization 
factor chosen here is different from that used in ^Sj. It does not effect any result but simply keeps 
the boundary weights finite in the limit Im(^) ±00. An obvious and important symmetry is 
given by ^ — — rA 



c± 1 
c± 1 



n, ^ 



B"-" c 



c±l 
c±l 



u, — r\ 



A special case is obtained setting a = 1 



B'^ ( r ± 1 



u, ^ 




r±l < «, 5 



1/2 

S/^i sm{^ lb u) sin(Ar + ^ =F u) 
sinA cosh2Im(0 



(2.5) 



(2.6) 



In this case we can use the equivalent form of (3.32) in JHI arising from its construction 

r-2 . 2 



r ± 1 



u, ^ 



sin A 



n 



sin"' A 



cosh2Im(^) J-jL sin(ti - {k + 1)A - ^) sm{u + kX + ^) 

r fl h /r-2 

r 

COS A 



m=2 



r-2 



\ m=l 



-(r-2)A 


-(r-3)A 






-(r-l)A 


-(r-2)A 




M-5-A 



(2.7) 



ei 



62 



er-2 
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As in |16j . we will vary the imaginary part Im(,^) of the field ^ to interpolate between the conformal 
fixed points. 

The given bulk and boundary weights satisfy the Yang-Baxter and Boundary Yang-Baxter 
equations and lead to commuting double row transfer matrices and an integrable model. The 
double row transfer matrices are also periodic, crossing symmetric and transpose symmetric 

Dna,Ha,(^,6,6) = Dnai|r,a,(^ + ^>ei,6), (2-8) 

Dna,Ha,(«,6,6) = Dnai|..a,(A-n,6,e2) (2.9) 

D.^a,Ha,(«,6,6) = {Ka,\r,a,iu,^l,^2)Y (2.10) 

for arbitrary complex values of the spectral and boundary parameters. 
Moreover, for ai = 02 = 1, the normalized transfer matrices 



t(u) = D.^^i|^,2i(u,^i,6)<S'r-i(tX,Cl)'S'r2(^^,6)'S'(^^) 

with 



sin(n -|- 2A) sin A "' 



sin(n -|- 3A) sin(ti -|- A) 



(2.11) 



^^""^ ~ sin(2tx + A)sin(2n-3A) ^ ^ 

Sr{u,^) = hr{u-Oh_riu + (2.13) 

sm A 

^^^^^ ^ sin A sin(n +{3- r)A) sin(n + {1 - r)A) ^4) 
sin u sm{u — A) sin(u -|- 2A) 



satisfy the functional equation 



t{u) t{u + \) = I + t{u + 3X) (2.15) 



where we have suppressed the ^ dependence. Periodicity, crossing symmetry and transpose sym- 
metry extend to this operator and its factors S{u), Sr{u,S,). The transfer matrix D is an entire 
function of u whereas the normalized transfer matrix t is a meromorphic function with poles aris- 
ing from the normalization factor. We need to determine appropriate choices of the boundaries to 
interpolate between distinct conformal boundary conditions of types (r, s) by varying the imagi- 
nary part of ^. We then need to solve the double row transfer matrix functional equation in the 
continuum scaling limit for the induced flow of the resulting eigenvalue (energy) spectra. 

In the scaling limit, the large corrections to the eigenvalues^ of the double row transfer 
matrix are related to the excitation energies of the associated perturbed conformal field theory by 



1 271" sm 1 

2log^nai|r2a2('") = ^ fb{u) + fr, a^lr^ aii^) + ai |r2 aa (0 + o(^ 



where /{, is the bulk free energy, friai\r2a2 is the surface (i.e. boundary dependent) free energy and 
Eriai\r2a2iO ^ Scaling function that, at the boundary critical points, reduces to 

Er,a,\r2a2L,^ = -^ + h + n, nGN (2.17) 

where h is one of the conformal weights allowed^ by the boundaries (?'i ai|?'2 02)- Observe that 
ai|r2 a2 '^2) in (|2.16|) Can be a function of up to two boundary parameters but we find that 
only one boundary parameter ^ needs to flow to reproduce the 931^3 thermal boundary flows. 

^We use D, t to indicate eigenvalues of D, t. 

^If the lattice is wrapped on a cylinder, its conformal partition function is given by a sum of characters determined 
by a "fusion" of the boundaries (ri ai\r2 02). 
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For our computations, we need certain information about the analyticity of the double row 
transfer matrices encoded in the zeros and poles of their eigenvalues. Because of integrability, 
D^^^l^^^^(ti, ^1,(^2) at different spectral parameters u,v forms a commuting family of operators so 
that its eigenstates can be chosen independent of the spectral parameter u and only the eigenvalues 
depend on u. These eigenvalues are Laurent polynomials in the variable z = exp(m). We have 
performed numerical diagonalization for small sizes (up to = 16 faces in a row). For these sizes, 
the zeros of the eigenvalues are extracted numerically and we extrapolate their pattern to the limit 
of large A^. We will refer to numerical observations obtained in this way as "numerics on D" . 



2.1 Critical point: classification of the zeros and TBA 



In the scaling limit at the isotropic point u = A/2, B^"" is critical if the boundary parameter is 
= A/2 or —A/2. These two choices yield the fixed and semi-fixed boundaries introduced in |21j . 
Actually, the parameter can be chosen in the regions ^ E [A/2, (5 — r — 1/2)A], ^ G [(1/2 — r)A, —A/2] 
respectively. Indeed, these are the maximal intervals where each single boundary contribution 
remains nonnegative at the isotropic point. The "numerics on D" shows that the scaling properties 
(in particular the classification of zeros that will be introduced shortly) do not depend on the actual 
value of ^ inside the intervals [211], so in this sense Re(^) is an irrelevant variable close to the critical 
point. 

In general, a cylinder partition function is a superposition of characters. However, a single 
character is obtained if the double row transfer matrix is built [501 with one r-type and one s-type 
boundary. The r-type boundary is obtained from B^^ by choosing ^ G [A/2, (5 — r — 1/2) A] or from 



^r+1,1 c];ioosing ^ E [—(1/2 + r)A, —A/2]. At the isotropic point u 
have explicitly (fixed boundary, "— " case) 



A/2 and with ^ = A/2 we 



B' 



rl 



B'-' r + 1 



A A 
2 ' 2 



A A 
2' 2 



r+1 




(2.18) 



sin \^JSr+l Sr 



Similarly, at ^ = —A/2 in B^^^'^ we have explicitly (fixed boundary, "+" case) 

r+1 

r + 1 ' 



j^r+1,1 



r + 1 



A _ A 
2 ' 2 




sinA-^<SV+i Sr 



^r+i.i lr + 2 



r + 1 
r + 1 



A _ A 
2 ' 2 



r+2 




(2.19) 



r+1 



In both the cases we are left with an alternating chain of r, r + 1 sites on the boundary. 

The s-type boundary weights are obtained as the braid limit of the r-type boundary weights 
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with r = s 

s 



lim B'^ s±l ^ 



s 



u, M = s±i- 



—2 sinA 

gi(2Re(?)+sA) 



(2.20) 



sin(^ + n) sin(A + — li) — 2sinA 



A similar expression holds for the limit Im(,^) — > +00. The u and Im(^) dependence in the last term 
cancels out explicitly. The overall factor is immaterial for our purposes because it corresponds to 
a trivial factor multiplying the transfer matrix. Among the r- and s-type boundaries there is the 
special common boundary condition with r = s = 1 corresponding to the vacuum. 

In the case of one r-type and one s-type boundary condition which leads to the single Virasoro 
character Xr,s{Q), a detailed analysis of the structure of zeros of the eigenvalues of the double row 
transfer matrices has been performed in [^. In the large N limit, each eigenvalue of the double 
row tranfer matrix is characterized by a specific pattern of zeros, organized as 1-strings or 2-strings 
in two analyticity strips PU] 

A 3A 

Strip 1: - 2 < < Y' ^*"P < ^^(""^ < '^^-^^^ 

The 1-strings are single zeros appearing in the center of each strip 

Re{u) = ^ or 3A (2.22) 

while the 2-strings are pairs of zeros {u,u') appearing on the edges of a strip, with the same 
imaginary part 

f (-A/2,3A/2), strip 1, 
(Re(u),Re(n'))= ' ' ^ . ^ (2.23) 

[ (2A,4A), strip 2. 

We use nil, "22 to denote the number of 1-strings in each strip and ni, n2 to denote the number of 
2-strings in each strip. For each (r, s) boundary condition these numbers satisfy particular (m, n) 
systems. The complete classification of the allowed patterns of 1- and 2-strings in terms of (m, n) 
systems and quantum numbers is summarized in Table [2 

For a generic function h{u), for example the normalized transfer matrix eigenvalues t{u), it 
is natural to introduce two functions /ii(x), h2{x) with reference to the center lines of the two 
analyticity strips 

hi{x) = + i^) for |Im(j;)| ^ vr (2.24) 

h2{x) = h{3X + i-) for |Im(3;)| ^ vr (2.25) 
5 

At the critical point of the TIM, the TBA equations derived in for the six (r, s) conformal 
boundary conditions in Table El are 

mi (1) _ 

ei{x) = - ^ log tanh(^L^) log(i + s2e-"2W) (2.26) 

k=l 

"12 (2) _ 

e2(x) = 4e-^-^logtanh(^L^)_i^*log(i + s^e-^i(^)) (2.27) 
k=l 



8 



Table 2. Classification, for all (r, s) boundary conditions of the TIM. of the allowed patterns of 1- and 
2-strings by (m,n) system and quantum numbers. The parity a ~ ±1 occurs when there are frozen zeros. 
The parities si,S2 = ±1 occur in the TBA equations. The expressions for ni, are only used on a finite 
lattice because in the scaling limit ni ^ N/2 — > oo. The number of faces in a row is even or odd according 
to N = {r- s) mod 2. 



Xr,s{Q) 


(m, n) system 


parities 


quantum numbers 


xi,i{q) 


mi, m2 even 

712 = mi/2 — 1712 

ni = [N + m2)/2 — rrii 


51 = 1 

52 = 1 


»,(-'^) n{ 7(1) _|_ ™, it'^ _l_ 1 m„ 

71^, — ^{J^j. + 777i — K) -\- I — 7712 

7lJ^^^ = 2(/f ^ + 7712 - A:) + 1 - 7711 




nil odd, 7712 even 

722 = {mi - a)/2 - 7712 

77,1 = (A^ + 7712 + <^)/2 — TTT-l 


51 = 1 

52 = -l 


0{ _|_ rn, i-^ _l_ 1 -m „ rr 

Tlj^ — ■^{J^j^ + TTll — K) -\- I — 7712 ~ 
7lS^^^ = 2(/f ^ + 7712 - A:) + 1 - 7711 + a 


X2,i{q) 


7771, 7772 odd 

772 = {mi + l)/2 - 7712 

77i = {N + 7772)/2 — 7711 


51 = -1 

52 = -1 


'){ tW _|_ ™_ _l_ 1 m„ 

71^, — ^V^fc + ?T7i — Kj + i — 7712 
7lJ^^^ = 2(/f ^ + 7712 - A:) + 1 - 7711 


X2,2{q) 


7711 even, 7772 odd 

772 = {mi - a + l)/2 - 7712 
77i = (A^ + 7n,2 + (T)/2 — 7771 


51 = -1 

52 = 1 


n^^"^ — 9(' r'-^-' 4- m, t"; 4- 1 m,, n- 
't^ — V fc ^ 1 y J- '"'2 — 'J 

nf'' = 2(/f ^ + 7772 - A:) + 1 - mi + a 


xz,i{q) 


mi ^ 2 even, 7712 odd 

712 = 77I1/2 — 7n,2 

m = (iV + 7772 + l)/2 - 7711 


51 = 1 

52 = -1 


n'j}^ = 2(4^^ + mi - A;) + 1 - m2 
71^^ = 2(/f ^ + 7712 - A") + 1 - mi 


X3,2{q) 


Till, 7772 odd 

772 = ("T-l - (T)/2 - 7712 

m = (iV + 7^2 + fJ + l)/2 - 7711 


51 = 1 

52 = 1 


= 2(lf ^ + mi - A;) + 1 - m2 - a 
77^^^^ = 2(/f ^ + m2 - A:) + 1 - mi + a 



where y^^ are the (scaled) locations of the 1-strings in the strip j = 1,2, that is to say, they are 
zeros of the normalized (and scaled) transfer matrices ti{x) = Sie~^^^^\ The parities are given 
in Table [21 The locations of the zeros are determined by a set of non- degenerate quantum numbers 
n) G Z by the quantization conditions 



= ''"''i'^' k = l,2,...,m2 (2.28) 



^2{y^''-i^) = iTrn'j!\ A: = 1,2,..., mi (2.29) 

The non-negative quantum numbers {ij^^} have the topological meaning that for a given 1-string 

(i) (i) (i) (i) 

y): , Ifl is the number of 2-strings with larger coordinate zf > y): ■ The conventional order of the 

zeros is y^^ < yg"*^ < . . . < ym]- The notation {I^\ ■ ■ ■ ,Iml\Ii^\ ■ ■ ■ ,Iml)a uniquely labels states. 
The corresponding expression for the scaling energy (|2.1H|) is given by j2()j 

= -Ye-y'^^^- ^e^Mog(l + S2e-^2(^)) (2.30) 



k=l 



c , (r — s){s — r — 1) m*Cm A m sr^sr^ ^(i) 
-24 + "- + ^ ^4 ^ + ^ ^ + 

j=i k=i 
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where C is the A2 Cartan matrix, A is a vector 



C 



2 -1 
-1 2 



1,-1 CT 



(2.31) 



and a = ±1 according to Tablel^Jif s = 2, 3 and cj = otherwise. The conformal partition function 
obtained from these energies yields a finitized fermionic character 



all states 



(2.32) 



24 1 '-rs I ■^ly^ • ^ 

mi,m2,o- 



mi + Til - 5a,l 
mi - 6„1 



m2 + n2 

7712 



(2.33) 



where the sum is over all states allowed by the finite (m, n) system and 

with the previous convention on a. 
2.2 Boundary flows 

The TIM admits 6 Cardy-type boundary conditions each corresponding to a single Virasoro charac- 
ter. Interpolating between them, as described in and shown in Figure ^ there are 7 integrable 
renormalization group flows, 5 generated by the perturbing operator 931^3 and 2 by (/Ji,2- In this 
paper we consider only those described by the A4 RSOS lattice model of ^H]) that is to say, the 
Lpi^3 flows. We expect that the ipi^2 boundary flows will be given by consideration of the dilute A 
RSOS lattice realization of the TIM. There is an important reason why A/i cannot describe both 
the perturbations, namely, the existence of two independent perturbing parameters would suggest 
the presence of an integrable surface instead of just the integrable line that is known to exist. 

From the A4 lattice description of the TIM we observe the presence of similarities between 
the ipi^3 flows related by inversion of the spin direction in the Affleck description. Moreover, 
from the conformal description, we also notice that three flows share the same "active" boundary, 
namely an s = 2 type boundary becomes an s = 1 type boundary, while the r-type boundary is 
just a spectator. We therefore divide the flows according to whether r remains constant or not 



constant r: 



Xi,2 = (0-) 

X3,2 = (0+) 
X2,2 = (d) 



Xi,i ^ (-) 

X3,l = (+) 
X2,l = (0) 



variable r: 



Xl,2 
X3,2 



(0-) 
(0+) 



X2A = (0) 
X2,l = (0) 



We will see that these similarities also arise in our TBA description of the flows. 



3 Boundary flows with constant r 

The common feature of this family of flows is the change from an s = 2 to an s = 1 label in the 
conformal character Xr,s{Q), keeping r fixed. Clearly, on the lattice this means that the boundary 
of type r is a spectator and the dynamics is given by the s type boundary. 

For convenience, from now on the complex boundary thermodynamic field on the lattice will 
be denoted ^latt- From (|2.2flj) we know that B^'^{u, ^latt) in the limit Im(^iatt) — ^ ice (and for 
all Re(^iatt) S M) is an s = 2 type boundary. We can take advantage of the fact that Re(^iatt) is 
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irrelevant at criticality to fix it to a convenient value. Using ()2.4|) and (|2.7p we have 



B^'M 2± 1 



att 



2±1 




2±1 




sin^ A 



sin(Ciatt + u) sin(A + ^latt 



' ' ' (3.1) 
and we see that if Re(^iatt) = — I™(6att) = 0, this boundary weight can be interpreted as the 
addition of a column to the lattice in presence of a B^'^ boundary that is of s = 1 type. So, on the 
lattice, the desired flow is driven by Im(^iatt) and can be summarized as 



S^'^(u,-A + iIm(6att)) 



Im(6att) 
Im(6att) 



iboo 
0, 



1 column added 



(3.2) 



We form a double row transfer matrix D^|2 i{u, —A + ilm(^iatt)) by coupling B^'^ with an r-type 
boundary on the left having no free boundary parameters. We are left with just one boundary 
parameter ^latt- 

By (|3.2|) . if Im(,^iatt) = ±oo, on the right we have an s = 2-type boundary coupled with an r 
type on the left, leading to the character Xr,2- If Ini(6att) = 0, the s = 1 boundary on the right is 
coupled with an r-type on the left, corresponding to Xr,i- We now introduce a simplified notation 
Ciatt = — A + with ^ real and 



|2,1 



[U, 



-A + 2^/5) 



(3.3) 



If we explicitly compute B'^'^{u, —A + i^/5) we see that it is real for arbitrary ^ and real u. The 
boundary weight on the left and the bulk weights are also real so each single entry of the transfer 
matrix is real for real u. This is nothing but the real analyticity property^ 



Dr:i|2,i(". -A + ie/5) = (D;^i|2,i(n*, -A + ie/5))*, for ah ^ 



(3.4) 



The same property holds for the normalized transfer matrix. Combining real analyticity, periodicity 
and crossing symmetry and using the notation introduced in (|2.24|) and (|2.25|1 we find the following 
reality condition for x G M 



^2 5^^ 

D(3A + 4,0 
5 



(D(^ + ^f,e))*, 

(D(3A + i|,e))*, 
tj{x,0*- 



(3.5) 

(3.6) 
(3.7) 



so that the transfer matrices are real along the center line of each strip. They also are transpose 
symmetric (|2.10|) : this implies that they are real symmetric and thus they have real eigenvalues 
leading to real energies. The previous equations are also true for the eigenvalues. Of course, if 
Re(.^iatt) / —A, the boundary weights are no longer real, the reality condition is lost, the transfer 
matrices become complex and so too do the scaling energies. Again, this is confirmed by the 
"numerics on D" . 

We now proceed with the flow xi.2 Xi,i which is the prototype among the flows in this 
family and so we describe it in detail. 



^We use * to indicate complex conjugation. 
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1 " 

1 



41 % 2X 3X 4\ 

B 

Fig. 2. The two mechanisms A, B respectively that change string content during the flow xi.i i-^- Xi,2- Note 
that this is the reverse of the physical flow. These mechanisms are illustrated in the upper-half w-plane for 
the states: A, ( ) (0)+; B, (0 0) ^ (0)_. 

3.1 RG flow xi,2 ^ 

From the previous discussion, the relevant double row transfer matrix is D(u, ^) = ^11^2 ~^ ~^ 
ii/'o) containing one boundary thermodynamic field ^. It is convenient to restrict to ^ ^ so that 
^ = —00 corresponds to the ultraviolet (UV) fixed point {xi,2) and ^ = to the infrared (IR) fixed 
point The adjacency rules force N to be odd. The addition of a column implemented by 

the boundary interaction (|3.1|) is consistent with the different parity of the number of faces in xi,2 
and Xi.ij as required by the adjacency rules 

N^^, odd 1 — > N^^ = N^^ + 1, even (3.8) 
3.1.1 Two mechanisms for changing the string content 

As in , it is sometimes convenient mathematically to consider flows which are the reverse of the 
physical flows. In this section the actual flow we consider is IR = xi,i ^ Xi,2 = UV. 

The pattern of zeros of the eigenvalues D{u, ^) change along the flow interpolating the extreme 
configurations described in Table |2 In particular, nii must change its parity from even to odd. 
Using (|2.2|) . a simple counting of powers of the variable exp(m) shows that D{u,(^) is a polynomial 
of order 2(A^ + 1) for ^ ^ (xi,i is included in this range) and becomes a polynomial of order 2N 
in the limit ^ — > —00. So, inside the periodicity strip —A < Re(n) ^ 4A, precisely two zeros must 
"escape to infinity" during the flow, one in the upper- and one in the lower-half u-plane. Indeed, 
from the "numerics on D" we find just two mechanisms for changing the string content during the 
flow. As shown in Figure |2 these involve the objects (1 or 2-string) which are furthest from the 
real axis in strip 1 at the IR point. Strip 2 is not directly involved. 

A. When the furthest object in strip 1 is a 2-string it moves away from the real axis, heading to 
+00; then it collapses into a pair of 1-strings; one of them disappears moving to +00 (it is 
actually reached only at the UV point), the other remains in the scaling region and becomes 
the top 1-string in strip 1. Of course, it cannot have a 2-string above it, so it is frozen and 
in the UV a = 1. This pair of 1-strings, called correlated 1-strings, has special features that 
will be explained in the numerical section. The disappearance of the top 2-string leads to a 
decrease by one unit in the quantum numbers in strip 1. 

B. If the furthest object in strip 1 is a 1-string it moves to -|-oo. The remaining 1-strings can be 
excited and are therefore not frozen, so in the UV we have a = —1. 

This qualitative description of the mechanisms implies the following changes in the values of the 
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quantum numbers 

1— > irii^ = + 1, a = +1 

A: » 4""' = 4"""-l. /l:;./™ = (3.9) 

^IR = n^^ - 1 

B: /;Jr,''' = 0: ^ ^yv ^ ^ir _ 1^ ^ ^ _1 (3 ^0) 

Observe that the change of parity of rrii is consistent with the known parities at the two endpoints of 
the flow, as in Table 121 The consistency and completeness of these two mechanisms is demonstrated 
by the explicit mapping between the finitized characters associated with the UV and IR fixed points, 
as shown in the next section. Moreover, we observe these same mechanisms when solving the TBA 
equations. 

The application of mechanisms A and B to the first few states is summarized in Table 01 



Table 3. Flow xi,i '-^ Xi,2 (reverse of the physical flow). We present the explicit mapping of states from IR 
to UV up to the UV level 6. Here n^^, n^V are the excitation levels above the ground states, respectively 
= and /i = 1/10. 



^IR 


Mapping of states - mechanism 




^IR 


Mapping of states - mechanism 







- (0)+ 


A 





6 


(40) ^ (4)„ 


B 


5 


2 


(0 0) ^ (0)_ 


B 


1 


6 


(31) ^ (2 0)+ 


A 


5 


3 


(10) ^ (1)_ 


B 


2 


7 


(1000|00) ^ (100|00)_ 


B 


5 


4 


(11) ^ (0 0)+ 


A 


3 


7 


(3 2) ^ (210)+ 


A 


6 


4 


(2 0) ^ (2)_ 


B 


3 


7 


(41) ^ (3 0)+ 


A 


6 


5 


(21) ^ (100)+ 


A 


4 


7 


(50) ^ (5)_ 


B 


6 


5 


(3 0) ^ (3)_ 


B 


4 


8 


(0 0) ^ (0 0)_ 


B 


6 


6 


(0000|00) ^ (000|00)_ 


B 


4 


8 


(1100|00) ^ (110|00)_ 


B 


6 


6 


(2 2) ^ (110)+ 


A 


5 


8 


(2000|00) ^ (200|00)_ 


B 


6 



3.1.2 RG mapping betvi^een finitized characters 

In this section we show that the two mechanisms A, B are compatible with the counting of states 
at the two endpoints of the flow, given by the finitized characters. 

We observe in ()3.9|) and 1)3. 10(1 that, starting in the IR, each state changes unambiguously by 
precisely one of the two possible mechanisms so that the counting of states is complete. More- 
over, the IR finitized character naturally splits into two terms corresponding precisely to the two 
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mechanisms A, B (see ()A.2|) ) 



(q) 



"12 +"2 

m2 



„IR 



m2+n2 

1712 



(3.11) 



24 \ g4"^ '^"^ 
B 



m2+ri2 
m2 



where we attach the label IR or UV only to the variables that change under the flow. The labels 
A and B on the sums indicate that the sums on m}^, m2 are restricted by the constraints imposed 
by the two mechanisms (|3.9|) and (|3.10j) . This can be understood by using (|A.3|) for the first strip 
in the second line of (|3.1H) so that the sum on rnf^ ^ 1 required by the mechanism A is manifest. 
Similarly, the factor related to the first strip in the last line of 1)3. 11() can be rewritten using (|A.4|) 
as a sum with the constraint Imi = required by mechanism B. 

The single IR energy level at the base of a tower of states fixed by the string content {mi, 1112) 
maps to a UV energy level according to 



A: 
B: 



q 24 qi^ '-'^ 



24 + 10 qA^ 



-i(myv_m2) 



24 qZ^ 



(3.12) 



(these mappings of energies are fixed by the known expression for the energies at the IR and UV 
critical points). Using 1)3. 9(1 and 1)3. 1U() . the mapping of the (7-binomials (counting polynomials) is 
^iven by 

-1 



A: 



B: 





1 












1 











Combining (|3.1H) to (|3.13)) we obtain the finitized UV character 



(3.13) 



Xi,i W 



24 ~ 10 



m2+n2 
m2 



(3.14) 



This shows that the completeness of the IR counting of states, together with the mechanisms A 
and B, implies the consistency and completeness of the UV counting. 



3.1.3 Solution of the functional equations 

To solve the functional equation for the eigenvalues of the double row transfer matrix D(ii, ^) = 
D^-^|2 ^{u, — A + i^/5), we follow steps similar to those in [201 but taking particular care in managing 
the terms containing ^. The functional equation can be rewritten using (|2.24|) and (|2.25j) as a 
coupled system between the two analyticity strips 

l + t2{x) (3.15) 
l + ti{x) (3.16) 



, .TT .TT 



^2(a; + ^-)^2(a:-^-) 
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1 1 

1 




1 1 
1 1 
1 1 

X O X / 
1 1 
1 1 
1 1 
1 1 
1 1 

xl/ ^ xl/ \ 




— 

1 
I 
1 

1 
I 
1 


— @ — I — > 

X 


2!^ 31 4^ 

1 1 
1 1 
1 1 

X O X 
1 1 
1 1 
1 1 
1 1 



Fig. 3. Analyticity, in the complex it-plane, of the order 1 contribution to the normalized transfer matrix 
eigenvalues t for the flow xi,2 '—^ Xi,i- Circles are zeros and crosses are poles. Double circles are order 2 
zeros. 



Im(rr)G(-^,^) (3.17) 

These equations can be solved by taking the Fourier transform of the logarithmic derivative of 
the equations, taking care to remove all of the zeros and poles that can generate a singularity in 
log tj{x). It is convenient to consider separately the order A^, order 1 and order 1/A^ contributions 
to tj{x). 

The analyticity of the order contribution is contained in the last factor of 1)2.11(1 , leading to 
a zero of order 2N at u = 3A and to two poles of the same order at n = 2A, 4A. These are at the 
edges of the second strip so they do not effect the RHS terms 1 + ti and 1 + 12 within (|3.17() . The 
order N contribution f{u) is thus given by the inversion relations 

/i(x + i|)/i(x-i|) = 1 (3.18) 
/2(x + i^)/2(x-i|) = 1 (3.19) 

The unique solution ^3] of these equations with the required analyticity properties is 

fi{x) = 1 (3.20) 
/2(x) = tanh^^J (3.21) 



If we divide (|3.15() . ()3.16|) by (|3.18() . (|3.19|) . we obtain a new system of equations containing the 
function ti{x) / fi{x) on the LHS and free of order zeros and poles in (|3.17|) 

-^(^ + i2);^(^-^2^ = l + t2{x) (3.22) 

f + = l + ti{x) (3.23) 

h 2/2 2 

The analyticity of the order 1 contribution in (|3.22|) . (|3.23|) . represented in Figure |S1 is specific 
to each boundary condition, being contained in the factors S'(n), Si{u) and S'2(ti, ^latt) of (|2.11|) . 
S{u) leads ti(x) to have a double zero at x = and two poles at x = itivr and leads t2{x) to have 
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a double zero at x = and two poles at x = itivr, that do not depend on the boundary parameter. 
Si{u) simplifies with an analogous term in the transfer matrix generated hy B^'^ yielding a constant 
term to the normalized transfer matrices tj{x). 5'2(u,^iatt) is crucial and specific for this family of 
boundary flows (in general it will be ^^(tt, ^latt) with Re(.^iatt) 7^ —A) 

art \ ^^"^^ sin(^iatt - u- X) sin(^iatt + u + 3X) cosh 2Im(^iatt) „ 

sm(4iatt - u) sm(^iatt - U-2X) sm(^iatt + u- X) sm(Ciatt +U + 2X) 

We see that S2{u,—X + i^/5) has zeros and poles in the second strip only: two single zeros at 
u = 3X± {x = and four single poles at « = 2A ± u = 4X± (x = ztivr it ^). The 
factor cosh(2Im(^iatt)) is compensated by the corresponding normalization assumed in l\2A\\ and 
always disappears from the equations. Observe that the given poles of order 1 are on the edges of 
the two analyticity strips so 1+^2(2;) and l + ti{x) in p.lSf) and H3.16() are free of poles inside ()3.17() . 
They are also non-zero in (|.S.17j) . To see this it is enough to observe that all factors tj / fj{x ^iir /2) 
have no order one zeros in |Im(a;)| < 7r/2. So we conclude that the order 1 contribution satisfies 

51(3; + i-^) 51(2; -i^) = 1 (3.25) 
92{x + i^) g2{x - i^) = 1 (3.26) 

The unique solution with the analyticity determined by the structure of zeros and poles just dis- 
cussed is 

5i(x,e) = tanh^l (3.27) 

52(2;, C) = tanh^ — tanh — ^ — tanh — ^ — (3.28) 

Clearly, the functions lj{u) = tj{u)/{fj{u)gj{u,£^)) giving the finite-size corrections satisfy the 
equations 

h{x + i-) h{x - i-) = l + t2{x) (3.29) 



2' ^ 2- 
.vr, 

'2- 



l2{x + i^)l2{x-i'^) = l + ti{x) (3.30) 



and do not possess the zeros and poles of order N and 1 that we have already removed. They 
only have other zeros, given by the 1-string zeros of D(ii, ^). Observe that these do not effect the 
functions 1 + t\{x) and 1 + t2{x), that again are non-zero on the real axis. Assuming N large, the 
position of the 1-string zeros v^^^ is real^ and we can remove their contribution by introducing the 
functions 

1 ■ 1 ■ 

Pji^) = n tanh[-(54^'^ - x)] tanh[-(54-''^ + x)], (3.31) 

k=l 



in each strip. It is simple to prove that these functions satisfy 

.vr .TT 



Pj{x + i^)pj{x-i^) = 1. (3.32) 



So, as before, we introduce the functions 



Observe that the classification of zeros in terms of 1- and 2-strings is only strictly true for A*" sufficiently large. 
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which are completely free of poles and zeros inside H3.17() . We can then solve the remaining equations 

h{x + i^)h{x-i^) = l + t2{x) (3.34) 

I2{x + i^)l2{x-i^) = l + h{x) (3.35) 

by taking the Fourier transform of the logarithmic derivative. The solution is 

log[i(x) = K*log{l + t2) + Cu (3.36) 

log[2(x) = i^*log(l + ti) + C2, (3.37) 

where Cj are integration constants. The kernel of the integration and the convolution are given by 

1 r+oo 

^(^) = n ^' {f*9){x)= dyf{x-y)g{y) (3.38) 

ivrcoshx J-oo 

Using (|3.33() we recombine all the contributions 

mi ^ 1 

logti(x) = log5ri(x,^) +^log [tanh[-(54^^ -x)]tanh[-(54^^ +x)]] 

k=l 

+K*log{l + t2) + Ci, (3.39) 

m2 



1 1 

^ogt2{x) = log/2(2;) + logg2{x,C) + ^log [tanh[-(57;f^ - x)] tanh[-(5vf ^ + x)]] 

k=l 

+A'*log(l + ti) + C2 (3.40) 



These are the lattice TBA equations. The boundary thermodynamic field enters the TBA equa- 
tions only through the boundary factors gi, g2- 

There are two important simplifications of the previous equations. It will be clear in the sequel 
that, in all the cases of interest, 

lim gj{x,^) = 1 or - 1. (3.41) 

X—^±00 

This expression and (|B.3|) from Appendix B can be used in the TBA equations (|3.39|) . 1)3.40(1 to 
compute the limit x — > +oo. This forces the integration constants to be an integer multiple of in 

Cj oc in (3.42) 

so that they can be removed from the equations and replaced by a sign, Sj = ±1, inside the LHS 
logarithm: log Sjtj{x). This amounts to fixing the branch of the logarithm. From (|3.7() we know 
that tj{x) is real for real x. From the analysis of the analyticity, we also know that 1 + tj{x) ^ 
for real x so we find that it never changes its sign. The limit in ()B.3|) forces it to be positive, 
1 + tj{x) > 0, so that the convolution term in the TBA equations is always real. 

3.1.4 Scaling limit and TBA equations 

The order 1/A^ term in (|2.16() . giving the scaling behaviour, is obtained in the scaling limit. Guided 
by the function /2, that has the well defined behaviour 

lim /2(x + logiV) = exp(-4e-^), (3.43) 

N—*oo 

we assume the general scaling forms 



hj{x)= lim hj{x + log N). (3.44) 
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This is consistent with the behaviour observed for the zeros of smah transfer matrices. Increasing 
N corresponds to adding 2-strings in the first strip, close to the real u-axis. This has the effect of 



pushing the original zeros v 



(i) ..,U) 



k ' 



W 



away with a logarithmic growth in their locations so that 



lim (5t;P^ - logiV), 



lim (5^^; 



(i) 



- v--fc -logiV) (3.45) 

are well defined. Applying these results to the lattice TBA system (|3.39|) . (|3.40|) we obtain the 
scaled-continuum TBA equations, after introducing the pseudoenergy functions 

ej {x) = - log Sjij (x) , i = 1, 2 

The properties of tj(x) discussed at the end of the previous section imply 



e-^i(^)gM, l + ij{x)>0, 
so that (for real x) the TBA equations are 

(1) 

X, 



ei[x 



e2ix) 



■mi (1) 

log giix,0-y, log(tanh 



k=l 



X G 



K * log |1 + S2e' 



(3.46) 
(3.47) 



4e - log52(x,0 - J]]log(tanh 



(2) 

Vk 



K * log |1 + sie 



(3.48) 



(3.49) 



k=l 



This means that any complex contribution to ej comes only from the source terms in the TBA equa- 
tions. Note that these perturbed TBA equations differ from the critical TBA equations ()2.26|) and 



p.27|) only through the appearance of the boundary perturbation terms log 51 (x, ^) and log 52 (2;, £,)■ 
The boundary thermodynamic field is also scaled so that it survives in the continuum limit 

5j(x, = lim gj{x + log iV, ^ - log N) (3.50) 



In the present case, we find from (|3.27|) . (|3.28|) that 

9i{x,i) = 1, 



tanh 



x + ^ 



(3.51) 
(3.52) 



so that only §2 remains in the equations. The sign used in the scaling of Qj^ for the boundary field, 
^ — log A^, corresponds to the value ^ = —00 at the UV conformal fixed point (xi,2(9)) and ^ = 00 
at the IR conformal fixed point (xi,i(9))- To be more general, we will keep track of gi in many of 
the equations that follow, because in the analysis of other flows this term will be non-trivial. 
In the TBA literature, it is standard to define 

Lj{x) = log |1 + ij{x)\ = log |1 Sje"'-'('')| G M. (3.53) 

The zeros of tj are then the zeros of Lj and vice versa, as is seen if we write the TBA equations in 
the more physical form 

mi 



Li{x) 



L2{x) 



log 



log 



1 -hsiffi(x,0 JJtanh 



k=l 



1 + e 



-Ae- 



's2g2{x,0 n ^^"^^ 



(2) 
Vk 



,K*L^{x) 



k=l 



(3.54) 



(3.55) 



Again following [20] , we can use (|3.3H|) , (|3.37j) to determine the scaling part of the eigenvalues 
of D(ti,,^) at the isotropic point u = A/2 



o »^i 



-Vk 



(1) 



A;=l 



pdx , , 

/ Z2^ ^2(X) 
J —00" 



(3.56) 
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3.1.5 Zeros and quantum numbers 



In the scahng hmit, the zeros of the double row transfer matrix satisfy 

ij{yl^^) = for the 1-strings (3.57) 
tj.(^(j) =p i7r) = for the 2-strings (3.58) 



The normahzed transfer matrix itself can have other order 1 zeros from the factors in (|2.1H) . In 
particular, the point x = —^^ being a zero for g2{x,S^) is also a zero for t2{x,(,), i.e. a zero in the 
center of the second strip. The string-like zeros can be obtained by a single equation; to see this, 
we perform the scaling limit and the shift x i— > x — i7r/2 in ()3.15|) . ()3.16|) 

ii{x) ii{x - iir) = l + i2{x-i^) (3.59) 

i2{x) i2{x - in) = l + ti(x-i|) (3.60) 
and find that each real solution of 

ij{x-i^) = -l (3.61) 

is either a 1-string or a 2-string in the strip 3 — j. The zero x = — ^ is not contained in the 
previous equation because it is coupled with a pole, as in Figure |2l The 1-string or 2-string cases 
are mutually exclusive, as the numerical analysis on small double rows transfer matrices confirm: 
the zeros generally repel one another and only single zeros are observed. Double zeros can appear, 
however, at exceptional values of ^, for example at the point where the 2-string in mechanism A 
transforms into the correlated 1-strings, as in Figure |2j 

The locations of the zeros are fixed by (|3.61|) so we can define the counting functions 

vr « vr 

ipjix) = —i ej(x — i—) = —i log Sjtj{x — i—) (3.62) 

that are a multiple of vr on a 1 or 2-string and introduce non-degenerate quantum numbers n^j^'^ G Z 
by the quantization conditions 

My'-k^) = -i £2(2/1'^ - i^) = n^k^TT, (^'^ + S2) = mod 2, (3.63) 
MVk^) = -i £1(2/^ - 4) = 4^^^' + si) = mod 2. (3.64) 



Note the inversion of the indices: ipi is for strip 2 and ip2 for strip 1. The same equations hold for 
the 2-string locations zf'^ so, each time ipj{x) = nvr is satisfied for an integer n with the appropriate 
parity, x is a 1-string or a 2-string in the strip 3—j. 
The explicit expressions for ipj are given by: 

V'i(x) = iloggi(x - i^,^) + i^logtanh( ^^ ^ 

k=l 



PV 



dy L2{y) 
27r sinh(x — y) ' 



(3.65) 



ma (2) 

'4)2{x) = 46""^ + i log 52(2; - + ^2_.^°Stanh(-^^ + i-, 

k=l 



- PV 



dy Li{y) 
2-K sinh(x — y) 



(3.66) 
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where the integral around the singularity x = y must be understood as a principal value. We take 
the fundamental branch for each logarithm so that in general log a + log b cannot be identified as 
log(a6). This prescription is convenient in view of the numerical computations and is used also in 
the TBA equations (HTTO . (HO^ . 

With this in mind, we can easily compute the limits of V'j for x ±00. For x —00, the 
dominant term in the expression of L2 and 1P2 is the exponential, so we have 

L2{x) ~ e"'^'^ — > 0, as X — > —00, (3.67) 

'il^2{x) ~ 4e~^ ^00, as 2; ^ —00 (3.68) 

These reflect the emergence of an infinite number of 2-strings in strip 1 as the size of the double 
row transfer matrix goes to infinity, as in Section 13.1.41 This dense filling of 2-strings forces the 
1-strings in the first strip only to be far apart from the —00 region, — co < y^^ < . . . < y'^^ so we 
can compute in full generality the limits 

lim Li(x) = log(l + si5i(-oo,0) = 1^°^^ (3.69) 

I— > — CXD I —00 

TT 

lim Vi(2;) = i\oggi{-oo-i-,C). (3.70) 

X— >— 00 Z 

In the present case, we have 

lim ^/;i(x)=0, for ^ > -00. (3.71) 

a;— >— 00 

The limit x +00 requires some care. We assume that Lj(+oo) is finite, as is confirmed by numer- 
ical computations, so that the contribution from the convolution disappears due to the principal 
value. In general, there are two possibilities 

{iloggj{+oo - if) - TTUij if yH] < +00, 

(3.72) 
iloggj{+oo - if) - 7r(mj - 1) if yli] = +00. 

In the present case, from the TBA computations at the IR critical point (xi,i('?))) we know that 
(7) 

y-mj < +00 in both strips so, away from the UV point, the previous equations reduce to 

tpj{+oo) = —TTrrij, ^ > —00. (3.73) 

Our main conclusion is that generally the functions tpj are globally decreasing, with nii, m2 
nonnegative and possibly large and again this is confirmed by numerics. There are a few exceptions 
that give rise to interesting behaviours that will be discussed later. The monotonicity is the 
ingredient we need to relate the two families {n^/^} and {ijf^} of non-degenerate and nonnegative 
quantum numbers. As in the critical case in Section 12.11 and in [2111, we start from the limit 
tpj{+oo). The 1 and 2-strings occur at multiples of vr with parity dictated by (|3.63|) and l|3.64() . so 
decreasing from one object to the next, the function -0^ increases by 2tt. The ordering of the zeros 
is < y2^ < . . . < ym] so above the zero y^^ there are /^"'^ 2-strings (by definition) and rrij — k 
1-strings. 

An expression for the parities Sj is still missing. If we assume the IR critical (xi,i(9)) TEA 
equations can be perturbed by the "small" function log ^2 = logtanh((x + ^)/2) with large we 
are led to maintain the even parity of mi, 772-2 and the parities si = S2 = 1, as in Table El This 
forces 

4^\4^^ odd (3.74) 
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We observed that ^/^j(+cxo) cannot be a zero. At least close to the IR region, we have y^"'^ < +00 in 
both strips, consistent with the limits 1)3. 73(1 where rrij are both even. The first available position 
for a zero is 'ipj{+oo) + vr = 7r(l — nij), yielding the following expressions for the quantum numbers 

= 2(4^^ +mi-k) + 1-7712, (3.75) 
raf^ = 2(4^)+m2-A;) + l-mi. (3.76) 

Observe that (|3.72j) is not predictive because we need to know in advance if a zero is at +00 to get 
the proper solution. 

A final comment on the zeros is useful. On the lattice, the 1- and 2-strings are zeros of the 
transfer matrix D{u,^). In addition to that, the normalized transfer matrix t{u,£,) can have other 
zeros and poles. After the scaling limit, only the functions tj{x,^) remain and, from the previous 
derivation, we see that t2{x,^) has a zero sX x = — ^. Restricting x to the real axis, the zeros 
of tj(x,£) and of Lj(x) coincide (|3.5.Sj) and can be either 1-strings or the boundary dependent 
zero X = — ^. In terms of the TBA functions Lj (or e^) and ipj, the characterization of the lattice 
quantities is: the 1-strings in strip j are the real zeros in common between Lj [x) and V's-j {x) — 27rn; 
the 2-strings in strip j are the real zeros of ip3-j{x) — 27rn and not of Lj{x); the boundary zeros 
are the real zeros of L2{x) and not of V'l — 27rn. 

As a complement to the description of the mechanisms and to the TBA equations, we sum- 
marize the main results concerning this flow starting from the IR point: 



mechanism A before the collapse and mechanism B: 

m 
2 



nil, fn2 even, n2 = ^ — m2 ^ 0, si = S2 = 1, 



n^^) = 2(/f ^ +mi-k) + l- m2, 
= 2(/f ^ + ma - A:) + 1 - mi; 

mechanism A after the collapse point: 

= nil + 2, = n^^\ = 1 — m2, 

the values of m2, n2, n^^"*, n^^^ with k < mf — 1, sj, remain unchanged. 

3.2 RG flow X3 2 ^ X3,i 

As pointed out in Section 12.21 this flow is dual to xi,2 > Xi,i ^-iid many of the considerations 
of Section l3.ll applv to it. For this flow we consider the double row transfer matrix Y){u,^) = 
■^3^i|2 i('")6att) with an odd number of faces in accord with the adjacency rules. We find the 
same change of parity in the number of columns ()3.8() and the same two mechanisms A and B 1)3. 9() , 
()3.1U|) changing the string content during the flow. The actual parities of N, nii are different but 
the changes are still fixed by the mechanisms A and B. The application of mechanisms A and B to 
the first few states is given in Tabled Again, we find a mapping between finitized characters 

^IR ^ ivUV = iviR_l 

showing the consistency of the mechanisms with the IR and UV counting as in Section 13.1.21 

In solving the functional equations, as noticed in Section r3.1.3l the order 1 behaviour is specific 
to each boundary condition. Here we need to consider S{u), Ss{u,X/2), 52(u,^iatt) from ()2.11() to 
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Table 4. Flow ^3^1 1-^ X3,2 (reverse of the physical flow). We present the explicit mapping of states from IR 
to UV up to the UV level 6. Here n^^, n^^ are the excitation levels above the ground states, respectively 
h = 3/2 and h = 3/5. 



^IR 


Mapping of states - mechanism 




^IR 


Mapping of states - mechanism 







[U U 




D 





5 


/A 1 

(41 


r\\ , ^ / Q n n 1 ^^ 
Uj 1-^ (oUU|Uj4. 


A 


5 


1 


(lU 




T5 

D 


1 


5 




Uj t-^ [Z L U|Uj-)_ 


A 

A 


5 


2 


/or* 

(2 U 


n\ /o 1 n\ 

Oj 1-^ (2 0)_ 


B 


2 


6 


/ 

(0000 


Ij 1-^ (000 lj_ 


B 


5 


2 


/'I 1 
(11 


Oj 1-^ (000 0j+ 


A 

A 


2 


6 


(10 


Oj 1-^ (100 0j_ 


B 


5 


3 


(3 




B 


3 


6 


(6 


0) ^ ((\\Q) 


B 


6 


3 


(21 


0) ^ (100|0)+ 


A 


3 


6 


(51 


0) ^ (400|0)+ 


A 


6 


4 


(40 


0) ^ (4|0)_ 


B 


4 


6 


(4 2 


0) ^ (310|0)+ 


A 


6 


4 


(31 


0) ^ (200|0)+ 


A 


4 


6 


(3 3 


0) ^ (220|0)+ 


A 


6 


4 


(2 2 


0) ^ (iio|o)+ 


A 


4 


7 


(1000 


1) ^ (100|1)_ 


B 


6 


5 


(0000 


0) ^ (000|0)_ 


B 


4 


7 


(2000 


0) ^ (200|0)_ 


B 


6 


5 


(5 


0) ^ (5|0)_ 


B 


5 


7 


(1100 


0) ^ (110|0)_ 


B 


6 



(|2.13|) . The only difference with Section ri.1.31 is the term 53(11, A/2) that replaces Si{u) 

Ssiu, A/2) = -_ ^ (3.78) 

sm(u + 2 j sm(n — ^A) 

It simply adds to Figure El poles at u = —A/2, u = 3A/2, that is, x = ^in. Being outside (|3.17j) 
they are not relevant, so we obtain the same order 1 system (|3.25j) . (|3.26|) . the same solution (|3.27j) . 
(|3.28|) and the same scaling terms p.5ip . 1)3. 52() : 

gi{x) = 1 (3.79) 
g2{x) = tanh^^^^ (3.80) 

The order 1/A^ analyticity is also specific to each boundary condition in the sense that the parity 
of rrij in (|3.3H) is boundary dependent, but the way to proceed is exactly the same as before. In 
conclusion, the equations obtained in Section [3. 1 .41 hold true also for this case. 
We summarize the description of this flow starting from the IR point: 

mechanism A before the collapse and mechanism B: 

nil ^ 2 even, m2 odd, ^^2 = ^ — "i2 ^0, si = 1, S2 = —1, 

= 2(/f ^ + mi - A;) + 1 - ma, 
nf) = 2(/f^+jn2-A;) + l-7ni; 



mechanism A after the collapse: 

mf^ = mi + 2, , = 7i^^\ = 1 — 1712, 

the values of m2, n2, with k < mf — 1, sj, remain unchanged. 
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Table 5. Flow X2,i X2,2 (reverse of the physical flow). We present the explicit mapping of states from IR 
to UV up to the UV level 6. Here n^^, n^^ are the excitation levels above the ground states, respectively 
h = 7/16 and h = 3/80. 



^IR 


Mapping of states - mechanism 




^IR 


Mapping 


I of states - mechanism 







(0|0) ^ (|0)_ 


B 





5 


(110|0) 


^ (11|0)_ 


B 


4 


1 




A 
A 


1 


5 


(100|1) 


1— > (i U| i j_ 


B 


4 


2 


(2\c]) (I oln), 

\ / \ / ~r 


A 


2 


5 


(5|0) 




A 


5 


3 


(000|0) ^ (00|0)_ 


B 


2 


6 


(300|0) 


^ (30|0)_ 


B 


5 


3 


(3|0) ^ (2 0|0)+ 


A 


3 


6 


(210 0) 


^ (21|0)_ 


B 


5 


4 


(100|0) ^ (10|0)_ 


B 


3 


6 


(200|1) 


^ (20|1)_ 


B 


5 


4 


(000|1) ^ (00|1)_ 


B 


3 


6 


(110|1) 


^ (11|1)- 


B 


5 


4 


(4|0) ^ (30|0)+ 


A 


4 


6 


(111|0) 


(0000|0)+ 


A 


5 


5 


(200|0) ^ (20|0)_ 


B 


4 


6 


(6|0) 


^ (5 0|0)+ 


A 


6 



3.3 RG flow X2,2 ^ X2,i 

The treatment of the flow X2,2 ^ X2,i is very similar to the previous flows. The relevant transfer 
matrix is now D(ii,^) = ^2i\2 i(^'^iatt) which has an even number of faces N. Again, we find the 
same change of parity in the number of columns ()3.8() and the same two mechanisms A and B (|3.9() , 
(|3.1())) . as illustrated in Table IH We also find a mapping between finitized characters 

xiT\l) ^ (3.81) 
^IR ^ ArUV^^iR_^ 

showing the consistency of the mechanisms with the IR and UV counting as in Section 13.1.21 

In solving the functional equations, we need to take care of the order 1 analyticity. With 
reference to Section [3.1.31 we need to consider the term 



S,{u, ^) = ;in^A^in(n + lA) sin(n + |A) 
^'2' sin2(u-iA)sin(n + fA)sin(u+|A) ^ ' 



All of its zeros and poles are independent of ^ and disappear in the scaling limit. Many of them 
disappear also at the lattice level: the factors sin(ti+|A) sin(n+|A) in the denominator are cancelled 
by equal terms in D(m,,^). Indeed, sin(u+ |A) explicitly appears in 52'i(ti, A/2); sin(n + |A) does 
not explicitly appear but we can prove its existence by crossing symmetry: D(u,,^) = D(A — u,^) 
implies that if u = — 3A/2 is a zero, then A— n = 5A/2 = — 5A/2 must also be a zero (the last equality 
is true because of periodicity). The remaining factors contain zeros at u = —A/2, u = 3A/2, which 
are irrelevant because they are outside of (|3.17|) . and a double pole at u = A/2 that compensates a 
double zero from S{u). So, already at the lattice level, gi{x) = 1. For §2, all the considerations in 
Section Hi. 1.31 can be repeated here. 

We conclude that the same TBA equations, quantization conditions and energy expression as 
in Section f3 . 1 .41 hold for this flow. 

We summarize the description of this case starting from the IR point: 
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mechanism A before the collapse and mechanism B: 

mi, m2 odd, n2 = "^"^^ — m2 ^ 0, si = S2 
4^) = 2(/f^+?ni-A;) + l-7n2, 
nf) = 2(/f)+m2-A;) + l-?ni; 

mechanism A after the collapse: 

= mi + 2, '^IIa.i = n^^A = 1 - m2. 



the values of m2, ^2, , with k < mf — 1, Sj, remain unchanged. 

4 Boundary flows with variable r 

4.1 RG flow X3,2 ^ X2,i 

For the flow X3,2 = Xi,3 ^ X2,i we consider the double row transfer matrix Y)^^^^{u,^\titt) with 
the trivial B^'^ boundary on the left, with no free parameter, and -B^'^(n, ^latt) on the right with 
the parameter ^latt £ C. The number of faces A'' is even. The limit Im(^iatt) —>■ ±oo on i?^'^(u, .^latt) 
reproduces the s = 3 boundary 1)2. 2Up and, for Im(,^iatt) = 0, i?'^'^(n, Re(^iatt)) yields an r = 2 type 
boundary (f2T8|) if Re(^iatt) G [-5/2A, -A/2]. If we choose ^att = -3A/2 + with ^ real. 



3 [S^ - sm(f A ± n + 4) sin(f A ± n - if) 

2 5y V ^3 sinA cosh2Im(0 



is real analytic in u. Notice that any choice other than Re(.^iatt) 7^ —3 A/2 would lead to the loss of 
real analyticity, that is a key property to get real flows. From the previous equation, we also see 
that each single entry of the transfer matrix is real for real u. This leads to the real analyticity for 
the whole matrix 

D(n,0 ^ D^[i|3_,(n, -^A + 4) = {B^^,\,^,{u* , -h + d))*, ^ G R. (4.2) 

The same property holds for the normalized transfer matrix. Combining real analyticity, periodicity 
and crossing symmetry and using the notation introduced in 1)2.241 I2.25() we obtain the reality 
conditions for x G M 

D(^ + z|,e) = (D(^ + .|,e))*, 

D(3A + i|,e) = (D(3A + i|,e))*, (4.3) 
t,(x,0 = tj{x,0\ J = 1,2. 

The transfer matrix is transpose symmetric ()2.1())) : with the previous equations, this implies that 
it is real symmetric so that its eigenvalues and the scaling energies are also real. The previous 
equations are also true for the corresponding eigenvalues. The integrable flow 



<i|3,(-,-^A + 4)-^ (4.4) 

X2,l if C = 



thus provides a lattice description of the renormalization group flow xi,3 ~* X2,i- The two ranges 
£ ^ and £ ^ actually describe the same flow. For later convenience, we choose £ ^ 0. 
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Table 6. Classification, for alternative realizations of (1, 3) and (2, 1) boundary conditions of the TIM not 
shown in Table |21 of the allowed patterns of 1- and 2-strings by (m,n) system and quantum numbers. The 
parity a = ±1 occurs when there are frozen zeros. The parities si,S2 = ±1 occur in the TBA equations. 
The expressions for ni, are only used on a finite lattice because in the scaling limit ni ~ N/2 oo. The 
number of faces in a row is even or odd according to N — {r ~ s) mod 2. The energy expression (j2.30|l still 
holds in these cases. 





(m, n) system 


parities 


quantum numbers 


N G even 


mi, m2 odd 

^2 = {mi - cr)/2 - 1712 

rii = {N + 1712 + — mi 


si = 1 

S2 = 1 


= 2(4^^ + mi - A:) + 1 - m2 - 
= 2(/f ^ + m2 - A:) + 1 - mi + a 


N G even 


mi odd, m2 even 

712 = (m-i + l)/2 - m2 

ni = (N + m2)/2 — mi 


51 = -1 

52 = 1 


= 2(4^^ + mi - A:) + 1 - m2 
nf^ = 2(/f ^ + m2 - A:) + 1 - mi 



If we compare the transfer matrix under consideration with those of Sections K-i.21 and we 
see that we are using inequivalent lattice realizations of the same conformal characters. Indeed, 
the parities and boundary spins of the corresponding chains are different 

Xl,3 = X3,2: D^^|3 -^(u,±icX)) / D^^|2 -^(^,±700) 

being even, A^' odd, the left transfer matrices are used in the present section and the right ones 
were used in the constant r flows. This means that we cannot use the rules given in Table |2J 

We used "numerics on D" to confirm the scaling behaviour in (|4.4|) and to classify the pattern 
of zeros appropriate for these lattice realizations. The results are shown in Table |21 We notice 
that for Xi,3{q) obtain the same {m,n) system as described in Tabled after a shift of the 
parity of the chain Ni^^ = N3^2 + 1- So, we have two different transfer matrices that give rise to 
the same finitization of the conformal character. In contrast, for X2,i(9) we obtain a new finitized 
character, built from a different (m, n) system. We argue that this must be the case from TBA 
considerations. Indeed, we expect two flows related by duality to enter X2,i- Because of duality, 
we expect the same boundary terms in the TBA equations, so the only way to generate different 
flows is to start from different and inequivalent patterns of zeros. 

4.1.1 Three mechanisms for changing the string content 

It is again convenient to consider the reversed (unphysical) flow, from IR = X2,i to UV = xi,3- 

The pattern of zeros of D{u,^) changes along the flow: mi remains odd while m2 changes 
from even to odd, according to Tabled D{u,^) is a polynomial of order 2N + 2 for > — oo in 
the variable exp(iu) inside the periodicity strip —A < Re(ti) ^ 4A and becomes of order 2N at the 
UV point itself. We therefore expect precisely two zeros to "escape to infinity" along the flow, one 
for each complex half plane (Im(n) < or lm{u) > 0). The "numerics on Z?" shows three active 
mechanisms that, for the most part, only involve the furthest zeros from the real axis in the second 
strip. They are distinct from the mechanisms described for the constant r flows but show many 
similarities. They were described in |16j for the flow xi,2 to X2,i, dual to the present case. More 
detailed "numerics on D" allows us to be more precise in describing the 2-string movements, as 
given in Figure E] and Table [3 
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Fig. 4. The three mechanisms A, B, C respectively that change string content during the flow X2.1 ^ X1.3 
which is the reverse of the physical flow. These mechanisms are illustrated for the states: A, (OOOjOO) i-^ 
(000|0)+; B, (111|00) ^ (000|1)_; C, (0) ^ (0|0)_. 



A. If the top object in both strips is a 1-string, the top 1-string in strip 2 moves to +cxo, decoupling 
from the system exactly at the UV point where it produces a frozen state with a = 1. 

B. If there is a top 2-string in strip 1 and a top 1-string in strip 2, they move away from the real 
axis towards +00. A 2-string comes in strip 2 becoming the top 2-string. Consequently, the 
quantum numbers in strip 1 decrease by 1 and those in strip 2 increase by 1. 

C. If there is a top 2-string in strip 2, it moves away from the real axis and collapses into a 
pair of correlated 1-strings. The furthest of these 1-strings from the real axis moves to -|-oo, 
reaching it precisely at the UV point. The 1-string closest to the real axis remains in the 
scaling region and becomes the top 1-string at the UV point. The quantum numbers in strip 2 
thus decrease by 1. 

These mechanisms are initially observed by the "numerics on D" then verified solving the TBA 
equations. For the mechanism B, the description of the movement of the 2-strings is still incomplete. 
We can imagine a first scenario where the relevant 2-strings remain exactly located on the border 
of the corresponding strip so that they need to reach -|-oo to be able to move from one strip to the 
other. If this is the case, a quite unusual exchange between the 1 and 2 strings in the second strip 
must occur. It can be avoided in a second scenario where the 1-string only reaches +00; the zeros 
forming the 2-string leave the first strip by the two sides pointing toward the second strip where 
they become the top 2-string. If this is the case, to keep track of the movement of the 2-string we 
need to solve a much more complicated system of functional equations involving 



ii{x)i2{x + i-tt) 



-TT < Im(x) < 0. 



A 2-string located between the first and second strip can be fixed by ^2(2; + ifvr) 



l + t2ix-i-) 



0, 



TT 



- < Im(x) < 



(4.6) 
leading to 

(4.7) 



where both the real and imaginary parts of x must be determined. This is equivalent to analytically 
continue ti, t2 out of the fundamental interval 1)3. 17() . Actually, the TBA equations are completely 
independent of the behaviour of the 2-strings so we can still compute scaling functions even if we 
miss the complete description of the movement of the 2-strings. 

The description of the mechanisms A,B,C forces the following changes for the various pa- 
rameters (if no label IR or UV is used the corresponding parameter does not change during the 
flow): 



A. /, 



(1) 

mi 



(2) 



IR 



0: 



IR 



m 



UV 



IR 



1, 



(4.. 
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B. /^\)''' > 0, I'^jJ'' = 0: 



m^2^ ^ m^^ = m^^ - 1, a = -1, 
/i'^ ^ Pk' -1, A; = l,...,mi, (4.9) 

/f" ^ /r" = /f" + l, ^ = l,...,m--l. 

C. /il'^ > 0: 

^IR ^ myv = miR + i, a = -l, /^^^^^^ = 0, 

if^" - '''' = - 1, fc=l,...,.n^ (4.10) 

The correlated 1-strings introduced here are in strip 2 while those introduced in Section fl-? . 1 . 1 1 were 
in strip 1 (and connected to the frozen case). 



Table 7. Flow X2,i ^ Xi,3 (reverse of the physical flow). We present the explicit mapping of states from IR 
to UV up to the UV level 6. Here n^^, are the excitation levels above the ground states, respectively 
h = 7/16 and h = 3/5. 



^IR 


Mapping 


of states - mechanism 


„uv 


^IR 


Mapping 


of states - mechanism 







(0) 


^ (0|0)_ 


C 





5 


(110|00) 


^ (110|0)+ 


A 


4 


1 


(1) 


- (1|0)- 


C 


1 


5 


(10 0) 


^ (100|0)_ 


C 


5 


2 


(2) 


^ (2|0)_ 


C 


2 


5 


(5) 


^ (5|0)_ 


C 


5 


3 


(000|00) 


^ (000|0)+ 


A 


2 


6 


(111|00) 


^ (000|1)_ 


B 


5 


3 


(3) 


^ (3|0)_ 


C 


3 


6 


(300|00) 


^ (300|0)+ 


A 


5 


4 


(100|00) 


^ (100|0)+ 


A 


3 


6 


(210|00) 


^ (210|0)+ 


A 


5 


4 


(0 0) 


^ (000|0)_ 


C 


4 


6 


(6) 


^ (6|0)„ 


C 


6 


4 


(4) 


^ (4|0)_ 


c 


4 


6 


(2 0) 


^ (200|0)_ 


C 


6 


5 


(200|00) 


^ (200|0)+ 


A 


4 


6 


(110) 


^ (110|0)_ 


C 


6 



4.1.2 RG mapping between finitized characters 

In this section we show that the three mechanisms A,B,C are compatible with the counting of 
states, given by the finitized characters [22], at the two conformal endpoints of the flow. 

From the definitions, we observe that for each IR state there is precisely one applicable mech- 
anism so that the counting of states is complete. Moreover, using ()A.2|) . the IR finitized character 
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naturally splits into three terms precisely associated with the three mechanisms 



mi , 









mi 







q 



mi-l+?i™' 



7ni — 1 



B 



jin^^Cm}^ mi 



mi+nY'-l 



+ E9 

B 



. A 



\m™CTn^^ mi 



'{N+ml^)/2-l' 




(mi-l)/2 


mi —1 







(Af+m^^)/2-l 



(W+mII^)/2 



(mi-l)/2 



(mi-l)/2 
^IR 



(4.11) 



Here and below we attach the labels IR, UV to the variables that change under the flow. The 
labels A,B,C on the sums indicate that the sums on mi, m}^ are restricted by the constraints 
imposed by the corresponding mechanism. The sum constrained by A can be understood using 
for both the g-binomial factors in (|4.1H) and corresponds to summing under the constraint 



I, 



(1) 



(2) 



IR 



m, — -I TR — 0. Similarly, the sum on B becomes apparent using (I A. 31) for the first strip q- 

binomial and ()A.4|) for the second strip g-binomial and likewise the sum on C uses (fO|l for the 
second strip only. The relations in Table IHl were used to rewrite the g-binomials as functions of 
mi, m^, N only. 

We now use the properties of the three mechanisms to map this expression for the finitized 
IR character (|4.11|) into the finitized UV character. An IR energy level at the base of a tower of 
states with string content fixed by (mi,m2^) maps to a UV energy level according to the energy 
expression (|2.3flj) that holds at the two conformal endpoints of the flow 



(4.12) 



^ ^Ai,3-igirr.UVcrr^UV^_l(„^_^yV) 
^ ^Ai,3-igirriUVCrriUV^+l(mi_rnUV)^, 
^A2,i - 5 glrnlRCmra^-m™ ^ ^Ai,3-i ^irriUVcmUV ^+ 1 („,j„„UV) ^ 

Using 1)4. 8() to (|4.1Up we rewrite the g-binomials appearing in the last three lines of 1)4.11(1 . Also 



A 
B 

C 
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taking into account the mapping of the energies H4.12() we obtain 



. A 







(mi-l)/2 


mi — 1 







(Af+mUV_i)/2 



(mi-l)/2 
mUV 



(mi-l)/2 



mi —1 



B 



irriUVCmUV (mi -mUV)^mUV 



1 uv " 




■™uv_^„uv_^- 


mi 




,nUV 



I UV 



UV T I uv 



(4.13) 



The last equahty in (|4.13|) follows from the qi-binomial identities in Appendix A. Again this argument 
shows the consistency of the IR and UV counting of states with the three mechanisms A,BiC. 



4.1.3 Order 1 analyticity and TBA equations 

The steps required to solve the functional equation for the scaling functions described in 

Sections 13. 1.31 to l3. 1.51 also apply for the flow xi,3 ^ X2,i-, the only differences being related to the 
boundary terms gj and the order 1 analyticity contained in the factors 5*1 (u), S{u), 5*3(7/, ^latt) of 
(|2.1H) . Si{u) cancels with terms from B^'^ and disappears. Because of the factor S{u), ti{x) and 
t2{x) both have a double zero at x = and two simple poles at rr = ztivr, independent of We 
now have 



'S'3(';^,6att) 



sinA sin(M + A + ^latt) sin(u + 3A - ^latt) cosh 2Im(^: 



lattj 



sin(n + ^latt) sin(u + 2A - ^att) sin(ti - A - ^latt) sin(u + 2A + ^latt) 



(4.14) 



where ^latt = — 3A/2 + so we find simple zeros at u = A/2 it that are in the first strip at 
X = ib,^, and simple poles at u = — A/2ibi^/5, 3A/2ibz^/5 that are in the first strip at x = zbivrit^ 
(all four combinations of signs). Remarkably, the second strip is free of zeros and poles dependent 
on ^. This leads to the same plot as in Figure 01 with the exchange of the two strips. Observe that 
the indicated poles are on the edges of the strips and the zeros in the middle so, using 1)3. 15(1 . 1)3. 16() 
we conclude that 1 + t\{x) and 1 + t2{x) are free of poles and are nonzero inside |Im(x)| < 7r/2 
()3.17|) . This implies that the order 1 factors can be removed from the functional system ()3.15|) . 
(I3.16|) using 



9i{x + i-)gi{x 



92{x + i-)g2{x 



I — 

2' 



(4.15) 
(4.16) 



with the indicated analytic properties. This system has the same form as for the constant r flows 
()3.25|) . 1)3. 26(1 . the difference being in the analytic content. The solution is the previous one (|3.27|) . 
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but with the exchange of the two strips 



51 (x,^) = tanh^ — tanh — ^ — tanh — ^ — (4-17) 
<72(x,e) = tanh^J (4.18) 



so the scahng hmit yields 



= tanh^ (4.19) 
52(^,0 = 1- (4.20) 

We are led to the same TBA equations ()3.48() . (|3.49() . energy expression 1)3. 56() and quantization 
conditions (|3.63|) to (|3.66)) as in the constant r flows. Similarly, the expressions for the quantum 
numbers (|3.75|) . (|3.76|) are not modified from the IR expressions so that 

mi odd, m2 even, si = —1, S2 = 1. (4-21) 

As a complement to the given description of the mechanisms and to the TBA equations, we 
summarize the main results concerning this flow starting from the IR point: 



mi odd, m2 even, n2 = ™V^"^ — m2 ^0, si = —1, S2 = 1, 



mechanism A and B, mechanism C before the collapse: 

'-11 

2 

2(/^'^+mi-A:) + l-7n2, 
^ = 2(/f ^ + ?n2 - A;) + 1 - mi; 

mechanism C after the collapse: 

A , o (2) (2) -I 

m^ = m2 + 2, n)Jj,_^ =n^jA = ^-mi, 



m. 



2 



the values of mi, ni, ra^ , with k < m2 — 1, Sj, remain unchanged. 
4.2 RG flow xi 2 ^ X2,i 

The preliminarly treatment of this case was presented in the letter |16j but we now have a more 
complete understanding of the mechanisms of movement of the zeros and of the role played by 
Re(^iatt) along the flow so we need to correct and clarify some statements. 

The analysis done in the previous section for the flow xi,3 X2,i can be repeated for this 
flow with D(n,.^) = ^ii\2 i(^j6att) and ^latt = I'^ + ^f - According to (|2.18|) and 1)2. 20() . this choice 
of the boundary field gives an r = 2 type boundary for ^ = and an s = 2 type boundary for 
^ — > ibcx). Moreover, this is the only value of Re(^iatt) that leads to real analytic transfer matrices 
()4.3|) and, consequently, to real positive energies. We need to analyse the order 1 analyticity. The 
relevant term is 

a ^ t \ s™-^ sin(^iatt -u- X) sin(^iatt + u + 3X) cosh 2Im(^iatt) / . „o^ 

sm(4iatt - u) sm(4iatt - U-2X) sm(4iatt + u- X) sm(4iatt + U + 2X) 

that has exacly the same zeros and poles obtained for the previous flow (|4.14j) : in the first strip 
there are two single zeros at n = A/2 it {x = it^) and four single poles at u = —A/2 it i^/5, 
u = 3A/2 lb i^/5 (x = zbivr it with no zeros or poles in the second strip. This leads to the same 
TBA equations and the same boundary term 

giix,0 = tanh^ (4.23) 
92{x,0 = 1 (4.24) 
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Table 8. Flow X2,i Xi,2 (reverse of the physical flow). We present the explicit mapping of states from IR 
to UV up to the UV level 6. Here n^^, n^^ are the excitation levels above the ground states, respectively 
h = 7/16 and h = 1/10. 



^IR 


Mapping of states - mechanism 




^IR 


Mapping 


I of states - mechanism 









A 





5 


(5|0) 


^ (4)_ 


B 


5 


1 


(1|0) ^ (0)„ 


B 


1 


5 


(110|0) 


^ (110)+ 


A 

A 


5 


2 


(2\0) (1) 


B 


2 


5 


(200|0) 


h->- (2 0)^ 


A 


5 


3 


(3|0) ^ (2)„ 


B 


3 


6 


(110|1) 


^ (110|00)_ 


C 


6 


3 


(0 0|0) ^ (0 0)+ 


A 


3 


6 


(200 1) 


^ (200|00)_ 


C 


6 


4 


(000|1) ^ (000|00)_ 


C 


4 


6 


(6|0) 


- (5)_ 


B 


6 


4 


(4|0) ^ (3)_ 


B 


4 


6 


(111|0) 


^ (0 0)_ 


B 


6 


4 


(100|0) ^ (100)+ 


A 


4 


6 


(210|0) 


^ (210)+ 


A 


6 


5 


(100|1) ^ (100|00)_ 


C 


5 


6 


(300|0) 


^ (3 0)+ 


A 


6 



obtained in Section 14.11 the only difference being in the following parameters that we take from 
|16j : mi, m2 are both odd, si = S2 = — 1 as at the IR starting point. 

The mechanisms were explained in |16j and the mapping between characters (see also Table EJ 

xZ\q)^xZ\'1) (4.25) 

was explicitly computed showing the consistency of the mechanisms with the IR and UV counting. 
Notice that tj{x) actually has no trace of the "pole" discussed erroneously there, instead the zero 
of 5i at X = — ^ gives a zero to ii(— and, correspondingly, a zero to Li(— ^). 

We summarize the main results concerning this flow starting from the IR point: 

mechanism A and B, mechanism C before the collapse: 

mi, m2 odd, n2 = ^^^^^ - m2 ^0, si = S2 = -1, 
n^^^ = 2(/f ^ + ?ni - A;) + 1 - jna, 
^ = 2(/f ^ + m2-k) + l-mi; 

mechanism C after the collapse: 

A , n (2) (2) , 

m, = m2 + 2, n 4 = n 4 = 1 — mi, 

^ ' mj— 1 jTij 

the values of mi, rii, n^^^ with k < ?ti^ — 1, sj, remain unchanged. 

5 Numerical solution of TBA equations 

The TBA equations (|3.48j) . (|3.49|) can be solved analytically at the two conformal fixed points, 
which occur at the endpoints ^ = ztco of the flow, leading to the results in [^. In particular, 
the energy expression (|2.3()|) and the finitized characters can be obtained analytically. It appears 
that the TBA equations cannot be solved analytically along the flow, however, with numerical 
computations we can obtain in detail the behaviour of the flow for intermediate values of ^ to 
interpolate between the known conformal fixed points. 
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Our numerical algorithm to solve the TBA equations and auxiliary equations for the zeros is 
iterative. We make suitable initial guesses for the pseudo-energies and the locations of the zeros 
close to one of the conformal fixed points. We then iteratively update in turn the pseudo-energies 
and locations of the zeros to find new values until the iteration scheme converges to the required 
accuracy. We then increment or decrement the value of the field ^ and repeat the process until we 
scan the full range in ^. This procedure is natural for the pseudo-energies in the TBA equations 
but iteration of the quantization conditions for the location of the zeros requires inversion of some 
complicated expressions associated with phases that wind along the flow. From the expression for 
ip2 l|3.63() . H3.66() . we can extract the zero y^^^ by one of the following: 

first strip algorithms: 

exp: inversion of the exponential term e~^; 

§2'. inversion of the boundary term, for the constant r flows only; 

1/^2 : scan the domain of the function 4'2{x) — T^n^P to find its zero; 

phase: inversion of one of the phases ilogtanh( t) used here). 

Similarly, from the expression for ipi H3.64() . (|3.65|) we can extract the zero by one of the 
following: 

second strip algorithms: 

gi: inversion of the boundary term, for the variable r flows only; 

ijji: scan the domain of the function 'tlJi{x) — -kii^ to find its zero; 

(2) 

phase: inversion of one of the phases ilogtanh(^^^ — ^). 

We start with the simplest choice, algorithm exp, phase in strips 1, 2 respectively. The computations 
at the critical points in [20] can all be done with this choice and do not require the introduction of 
other algorithms. We emphasize that these different algorithms are not equivalent under iteration 
and it can happen that some schemes fail to converge in some intervals. In every case, however, 
the validity of our use of different algorithms in different intervals is confirmed by the fact that 
the plots of the scaling energies, and the locations of all of their associated zeros, join smoothly in 
some typically small overlapping interval of convergence. 

The contraction mapping theorem, says that if a mapping / : M — > M on a complete metric 
space M is a contraction then there exists a unique fixed point xq = /(xq) and all the sequences 
obtained under iteration starting from an arbitrary initial point x £ Ad converge to the fixed point. 
For our purposes, we take M C R to be a closed interval and a contraction means that the first 
derivative is ^ A; < 1 inside M, for a fixed k. In particular, |/'(xo)| is a measure of the rate 

of convergence. We will see that, in some cases, [/'(xq)] ^ 1 and although the fixed point exists it 
cannot be reached by iterating the mapping /. 

5.1 Flow xi,2 ^ Xi,i 

The scaling energies for the states in Table |31 are presented in Figure El The second strip is a 
spectator and is not effected by mechanisms A and B in 1)3. 9|) . 1)3. 1U() . so it always has m2 even and 
is not relevant in the following discussion. As usual, for numerical convenience, we describe the 
reverse of the physical flow. 

From Table |31 we begin by following the simpler mechanism B that flrst occurs in the 
excitation (0 0) (0)_ and has only one strip of zeros. We start decreasing ^ from the IR point 
^ +00 with algorithm exp and make the following observations: 

• the numerical iteration converges in the interval ^ ^ — 2 but fails for smaller values; 



32 



-10 



-5 







5 



Fig. 5. Scaling energies for the flow xi,2 '-^ Xi,i- The list of states is given in Table El The intermediate 
region of the mechanism A levels (shown dashed) are schematic and have not been obtained from the solution 
of the TBA equations. 
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• the zero 1/2 moves towards +00, consistent with the lattice description of mechanism B while 
moves slowly in the finite region, as shown in Figure El 

If we denote by f{x) the mapping obtained by inverting the exponential term in (|3.66)) to solve 
for X and follow the value of f'{y^2^) as we vary ^ we see that f'{y!^^) < 0.1 for large ^ > but 
f'{y^^) —00 for ^ — > — 00 while it remains bounded for the other zero, ^0.1, for all 

values of ^. This is a clear indication that the zero y^^ = ym\ cannot be reached by iteration. We 
certainly know that it exists because in the whole interval ^ > — 2 we can track it and at ^ = — 2 
it is at a finite position, y^^ ~ 1.49. An analytic estimation of f'{y^^) shows that the dominant 

behaviour is f'{y2^) ~ — e~^2 ' if y^'^ is sufficiently large, as we expect in the region ^ — > — cx), 
see also Figure IHl These considerations show that we need a different algorithm to determine the 
location of this zero. 

If we try with algorithm 52 we observe the opposite behaviour for the mapping f{x) obtained 
by inverting the boundary term log<^2 in l!3.66j) : f'{y^^) —>■ —00 for large ^ » and f'{y^^) is 
small for ^ < —1.5 so we can start from the UV point and combine the results from this algorithm 
with the previous ones, thanks to the existence of a small overlapping region of applicability. 

Actually, we find that the algorithm -02 is more general in the sense that it can be applied 
everywhere, being intrinsically non- iterative. We always use this choice for the largest first strip 
zero ym| (the other first strip zeros can be fixed with the slightly faster algorithm exp). A similar 

problem shows up even in the second strip for higher excited states so we use the non-iterative 

(2) 

algorithm ipi for all the second strip zeros y^, , = 1, . . . , m2. 

Once the appropriate algorithms are applied, the numerical analysis of the states proceeds 
smoothly from IR to UV and the predictions of mechanism B are completely confirmed. In par- 
ticular, the IR description of zeros and quantum numbers can be used all along the flow and ymj 
only reaches +00 at the UV critical point. 



6 - 



4 - 



2 - 



- 




-2 I ' ' ' ' 1 

-10 -5 5 10 

Fig. 6. Movement of the zeros versus ^ for the state (0 0|0 0)_ 1-^ (0 0|0 0), described by mechanism B, 
in the flow xi,2 ^ Xi,i- Solid lines from bottom to top are . . . , 1/4^'', dashed lines from bottom to top 
are y^^' . This picture is consistent with the lattice predictions in part B of Fig.|5| 
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zeros 




I ' ^ 

-20 -10 10 

Fig. 7. Scaling energy and movement of the zeros for the state (0)+ ^ { ), involving mechanism A, in the 
flow xi,2 '—^ Xi,i- lower plot, the solid line corresponds to the position of the correlated 1-strings and 

the dashed one to the top 2-string in strip 1. 



The mechanism A requires some care. We examine the ground state ( ), as in Table |21 
starting in the IR where there are no zeros and decrease This procedure must fail somewhere 
because we expect the top 2-string in strip 1 to transform into a pair of correlated 1-strings, from 
Figure l2j Keeping track of the top 2-string (with algorithm exp) we see that the TBA equations 
reproduce exactly the movement toward infinity, as shown in Figure [7| up to a value ^ > —6.5. 
We next start from the other end of the flow, the UV point, now increasing ^ from — oo. There 
and only there, one of the correlated 1-strings sits exactly at +oo and is essentially decoupled from 
the system, restoring the correct odd parity of mi. Away from the UV point we need to consider 
both of the correlated 1-strings, so now mi = 2. To get the code working, we need to use the 
algorithm g2 for the largest zero (?/mi in general) and use the usual algorithm exp for the other 
zero (k = 1, . . . ,mi — 1, in general). We reproduce the expected behaviour, as shown in Figure [7| 
for e < -9.5. 

In the interval ^ G (—9.5, —6.5) containing the collapse region the TBA equations are difficult 
to solve numerically. Indeed, the indicated zeros are not in the standard positions; in terms of 
u their real part is not fixed and must be determined as well as the imaginary part. In solving 
(|c{.63|l to obtain the position of the two 1-strings for ^ < —9.5 we see they have the same quantum 
number = n^2^ = 1 as they had when they were a 2-string (^ > —6.5). This explains the name 
correlated 1-strings because they are on a different footing from the ordinary 1 and 2-strings. At 
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criticality, the furthest among them is at +00 and the closest is the so called frozen zero already 
observed in [2111 ■ The monotonic decreasing behaviour of the function ip2ix) assumed to introduce 
the quantum numbers is observed to fail in the present case, in a neighbourhood of the furthest 
zero among the correlated 1-strings. In this way the space to allocate two degenerate quantum 
numbers is created. We see the same behaviour in all the mechanism A states. 

We have a consistency check. At the critical point Xi,2, using the data in Table 121 we have mi 
odd, 777-2 even, even and n^^^ odd. Comparing with the the limits 1)3.72(1 . we have agreement 
only if ^pj(-'t-oo) = —Trrrij and ym] < +00. For the first strip, 'ip2{+oo) = —m2TT is an even multiple 

of vr and n^^^ is even so, in addition to the known zeros labelled 1, . . . , mi, we have room for another 
zero located at +00. This zero was not counted in [20] because it is "almost" decoupled from the 
system but now we know exactly its origin and its role from the off critical behaviour. In the 
mechanism A case, that is the frozen case {a = +1), this zero is the partner of the frozen one in 
the correlated 1-strings; in the mechanism B case, corresponding to the unfrozen case (fi = — 1), it 
is the zero that left the finite region; in this case it is a normal 1-string. 

Repeating the previous argument for the second strip, it seems that we have room for an 
additional zero. Actually, it is only there at the conformal point xi,2- during the flow it cannot 
exist either at +00 or at a finite location so we don't need to consider its presence at all. 

5.2 Flows X3,2 ^ X3,i and X2,2 ^ X2,i 

These flows behave similarly to xi,2 Xi,i described in the previous section and the algorithms 

to solve the TBA system are exactly the same. In Figures |H1 and El we plot the scaling energies for 

the states in Tables HI and respectively. They are complete up to the top level in each plot. 

A comment for the second of these flows is useful. At the IR point and all along the flow except 

(2) 

at the UV point, the value of is even and, in addition, 7/^1 (—00) = is even ()3.71|) so that at 

(2) 

X = —00 there can be a zero y\ . Actually, a 1-string occurs whenever its quantum number equals 

(2) 

the number of 2-strings in the same strip, = 77-2. In this case, the zero sticks to its position for 
all values of and the TBA equations can be slightly simplified taking the appropriate limit. 

5.3 Flows xi,3 ^ X2,i and xi,2 ^ X2,i 

We consider first xi,3 ^ X2,i, displayed in Figure El It doesn't present any new difficulties with 
respect to the description in Section [5.11 In this case, we have three algorithms for each strip to 
get the position of the zeros and we need to choose those that converge in the appropriate regimes 
and depending on the mechanism. 

In general the first strip zeros labelled by 1, . . . ,7771 — 1 do not require special care and the 
simplest algorithm exp can be applied to them. So, in the following discussion, we refer just to the 
remaining first strip zero 7771 and the second strip zeros. In general, they can be determined by the 
algorithms ipj, j = 1,2 except for a few exceptions. 

Considering mechanism A, we just have a zero escaping to +00 which causes no problems 
for the numerical solution. The fiow of these states can be followed without interruptions from IR 
to UV. 

Considering mechanism B, we start from the IR and observe that the movement of the 2- 

strings depicted in Figure 21 doesn't influence the TBA equations, as expected because they are 

just spectators; in fact the scaling energies are determined by the 1-string dynamics only. The zero 

(2) 

moving toward infinity, y It,, is most conveniently determined by inverting the boundary term gi. 

Considering mechanism C, we start from the IR and we expect the iteration to fail because 
new objects must appear in the region where the transformation of the zeros depicted in Figure jlp 
occurs, namely —9.5 ^ ^ —6.5. Data for smaller values < —9.5 can be obtained by starting 
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I \ \ \ 1 1 

-12 -8 -4 4 8 

Fig. 8. Scaling energies for the flow X3,2 '-^ X3,i- The list of states is given in Tabled The intermediate 
region of the mechanism A levels (shown dashed) are schematic and have not been obtained from the solution 
of the TBA equations. 
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-12 -8 -4 4 8 

Fig. 9. Scaling energies for the flow X2,2 '-^ X2,i- The list of states is given in Table [3 The intermediate 
region of the mechanism A levels (shown dashed) are schematic and have not been obtained from the solution 
of the TBA equations. 



38 



I 1 \ \ 1 1 

-10 -5 5 10 

Fig. 10. Scaling energies for the flow xi^3 i-^ X2,i- The hst of states is given in Table [7| The intermediate 
region of the mechanism C levels (shown dashed) are schematic and have not been obtained from the solution 
of the TBA equations. 
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Fig. 11. Scaling energies for the flow xi,2 '-^ X2,i- The hst of states is given in Table|Sl The intermediate 
region of the mechanism C levels (shown dashed) are schematic and have not been obtained from the solution 
of the TBA equations. 
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from the UV point and including the zero at infinity, y = +00. Again, this particular zero is 

most conveniently determined by inverting the boundary term gi. 

A preliminary discussion of the flow xi,2 ^ X2,i was given in JB]. It is very similar to the flow 
Xi,3 X2,i- For completeness, we present our current, more accurate numerical results in Fig. 11. 



6 Conclusions 

In this paper we have used a lattice approach to derive exact TBA equations for all excitations of 
the 5 (fi^s integrable boundary flows of the tricritical Ising model. We have shown that, along these 
boundary flows, the patterns of zeros classifying the states can change by one of 2 or 3 mechanisms 
which have been explicitly identified. These mechanisms produce precise mappings between the 
relevant finitized characters describing the finitized energy spectra at the conformal UV and IR fixed 
points. The TBA equations were also solved numerically to determine the interpolating boundary 
flows for the leading excitations. 

Even for the TIM there remain other questions of interest. First, it is desirable to have a 
comparison of our results with the results of the TCSA. The methods used here could usefully 
be generalized to study other flows, such as X2,2(9) Xi,ii^) ~^ X3,i(9)) which involve linear 
combinations of Virasoro characters. It would also be of interest to study the missing 2 ipi^2 
integrable boundary flows, but presumably, this would involve a study of the relevant dilute A 
lattice model. Likewise, our approach should be applied to the boundary flows of the TIM from the 
superconformal perspective, especially since the integrable boundary conditions on the lattice that 
correspond to the superconformal boundary conditions in the continuum scaling limit are already 
known j23j . Similarly, it remains to incorporate the flow of boundary entropies into our approach. 
Ultimately, of course, the considerations of this paper should be generalized to all the minimal 
models. 
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Appendix A. Gaussian polynomials 

The expressions for the finitized characters are based on the following relations for the (7-binomials: 

n 1 1 Im-l 
/l =0/2=0 lm=0 



Lm J 



rn — 1 



(A.2) 



The following expressions are useful to understand the contribution of each mechanism to the full 
character. Indeed, the restriction > leads to 



n I\ Im-l 
h = lh = l Im = l 



Jl 



. + /„ 



n-1 I[ 

/( =0/^=0 i'm=o 



m+n— 1 



(A.3) 
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whereas the restriction /„ 



= leads to 

n II Im-2 
/l=0/2=0 Im-1=0 



,/l + ... + /„i_l 



m— 1+n 
m — 1 



(A.4) 



Appendix B. Braid limit 

In this section we consider the braid hmit of the normahzed transfer matrix and its eigenvalues. 
By direct computation, we obtain the following limits: 



lim e 

Im(M)— >±oo 











X-u, g / 


9 




u 



lim 

Im(it)— >±oo 




4sin^A 



(B.l) 



lim e 

Im(ii)— >±oo 



±4iu 




,±2iA 



cos A 



8 sin"* A sin sin 2^2 



(B.2) 



We can recursively apply the first result to the normalized transfer matrix to remove all the columns 
then use the second expression to take care of the boundary terms. This leads to the conclusion 
that t(u,^i,^2) is diagonal and the actual value of the limit is the same for all eigenvalues and all 
integrable boundary conditions, that is to say that it is proportional to the identity. For a generic 
eigenvalue we thus have 



lim t{u, ^,1,^,2) = 2 cos A 

Im(«)— ►itoo 



1 + V5 



(B.3) 



In j20j the same expression was boundary dependent because of different normalizations in the 
right boundary term. 
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